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Shuang Xu and Yang Yvfl 

National Laboratory of Solid State Micro structures and Department of Physics, Nanjing University, Nanjing 210093, China 



o 

(N 
Oh' 

(N 
(N 

S 1 

9 L,: 
-i— > 

cr 



> 
in 

On 

ON 

O 

o 



X 



Guozhu Sun 

Institute of Superconductor Electronics and Department of Electronic 
Science and Engineering, Nanjing University, Nanjing 210093, China 
Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA 

(Dated: September 23, 2010) 

We proposed a new method to measure the energy spectrum of a superconducting flux qubit. 
Different from the conventional frequency spectroscopy, a short triangle pulse is used to drive the 
qubit through the anticrossing and generates Landau-Zener-Stiickelberg interference patterns, from 
which the information of the energy spectrum can be extracted. Without installing microwave 
lines one can simplify the experimental setup and reduce the unwanted effects of noise. Moreover, 
the method can be applied to other quantum systems, opening the possibility of calibrating and 
manipulating qubits with linear pulses. 



PACS numbers : 85.25.Cp, 03.67.Lx 



INTRODUCTION 



It has been shown, both theoretically and experimen- 
tally, that superconducting devices based on Joseph- 
son junctions are ideal architectures successfully display- 
ing macroscopic quantum phenomenal. On Josephson 
junction devices, typical quantum phenomena like quan- 
tum tunneling 2 - - — , energy level quantizations, and coher- 
ent superpositions^— have been predicted and experi- 
mentally demonstrated. Temporal evolution of quantum 
states and Rabi oscillations have been observed in macro- 
scopic sens o 10 i n . By virtue of the convenience in engi- 
neering, many new ideas are expected to be implemented 
with Josephson junction devices; one of them is quan- 
tum computation^^. The basic requirement of quan- 
tum computation is a quantum bit (qubit). In practice, 
accurate, efficient and stable manipulations on qubits 
are the prerequisites of quantum computation. People 
have proved Josephson Junction devices as a promising 
candidate for a qubft^ i 11 ' 13 ! 14 . Furthermore, integrated 
multi-qubit control has been attempted with Josephson 
Junction devices 1 ^— . 

In order to control a qubit with high fidelity, first 
of all, we have to measure the energy spectrum of 
the qubit. Conventional spectroscopy utilizes frequency 
resonance to measure the absorption spectrum of the 
qubit. This method originates from the historical fre- 
quency spectroscopy that has been applied to atoms or 
molecule a 20 i 21 . In the case of Josephson junction devices, 
which are often referred as artificial "atoms" , people in- 
ject microwaves to the junction circuit and measure the 
absorption spectrum. However, this method becomes 
challenging if the energy is over 10 GHz, because mi- 
crowaves of this high frequency are hard to optimize with 
common electronic devices. Although high frequency mi- 
crowave sources are available, they are expensive and the 
signal becomes noisy after passing multipliers. Transmis- 
sion lines and waveguides also limit the applications due 



to restricted bandwidth 2 ^. To overcome this challenge, 
some groups have developed alternative techniques such 
as photon-assisted tunneling 2 ^. In 2008, Berns et ai. pro- 
posed a method called amplitude spectroscopy^. They 
drove a flux qubit with a microwave of 0.16 GHz, whose 
amplitude was as large as to sweep through anticross- 
ings in higher energy levels and produced Landau-Zener- 
Stiickelberg (LZS) interference. They observed diamond- 
like interference patterns, from which the information of 
energy levels with energies up to 100 GHz was extracted. 
However, microwaves have phase indeterminacy, which 
introduces decoherence of the qubit and leads to diffi- 
culties in exact manipulations. So we propose another 
method replacing microwaves with a short triangle pulse. 
We have studies the LZS interference generated by a tri- 
angle pulse sweep. With simulations we have obtained 
the LZS interference patterns, which encode the informa- 
tion of the spectrum. Remarkably, this technique is also 
useful in coherent manipulations on quantum systems be- 
cause the signal length can be set within a sufficiently 
short time. 

This article is organized as follows. In Section |H] we 
shall introduce the model in our simulations. Section IHlI 
is devoted to the results and analytical explanations of a 
two- level system. In Section HVl we take a further step to 
study the case of a multi-level system. Conclusions are 
made in Section fVl 



II. MODEL 

The qubit under our discussion is similar to the one in 
the experiment by Berns et al£±. It is a superconducting 
flux qubit built with a superconducting loop interrupted 
by three junctions. The qubit is biased with an exter- 
nal flux $ ex t around <E>o/2, where $o = h/2e is the flux 
quantum. The potential of the qubit has the double-well 
shape parameterized by the flux bias. According to quan- 



2 



turn mechanics, there are quantized energy levels located 
in both wells. If the noises are suppressed to a sufficiently 
low degree (much smaller than the energy gaps), quan- 
tum behaviors can be observed. We label the localized 
states as \Ln) (left well) and \Rn) (right well), where 
ri = 0, 1, 2..., counting from the ground state. In the ba- 
sis composed of the localized states, the Hamiltonian of 
the qubit has non-zero off-diagonal elements, which mark 
the interwell tunneling. Diagonalizing the Hamiltonian 
and calculating its eigenvalues, we obtain the spectrum 
of the qubit theoretically. The spectrum is often param- 
eterized by the flux detuning Q ex t — ^o/2, as drawn in 
Fig-Q]). There are anticrossings at certain values of the 
flux bias, avoiding energy degeneracy of the two states in 
different wells when they arc tilted to the same energy by 
the external flux bias. In regions other than the vicini- 
ties of the anticrossings, the energy varies linearly with 
the flux bias, which is a valid approximation in common 
cases. 
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FIG. 1: (color online) Typical energy-level diagram of a flux 
qubit, showing the relation between the energy and the flux 
detuning <& ext - $ /2. 

The study of Landau-Zener transitions and 
Stiickelberg-typc interference in a qubit has lasted 
for a long time^ir— . Usually, the experiments are 
carried out as follows. The qubit is initialized at the 
ground state in a well, e.g. |_L0) . Then a microwave 
signal is injected through the bias line. After some 
periods of sweeping, the population in a certain state 
is measured by using a dc SQUID juxtaposed aside 
the qubit to detect the loop current direction, which 
is decided by the state of the qubit. The previous 
works have observed clear interference patterns and 
managed to extract useful data from the patterns. In 
our simulation, we follow the same key steps, but have 
changed some details. First, our sweep signal is a linear 
triangle pulse, which can be accurately programmed 
and generated. This seems a minor change but brings 
great convenience to the experimenters. Moreover, this 
is closer to the type of transitions originally studied by 
Zener so Zener's formula^ is of higher accuracy in our 
work. Second, we measure the system every time when 
a single pulse ends. This ensures that we will obtain 



the direct result of one-turn interference within the 
decoherence time of the system. The result has not been 
disturbed or diminished by effects of the environment 
or defects. Therefore, it is a convenient operation with 
both accuracy and efficiency. 

Now we simulate the revised process of the LZS in- 
terference. We drive the qubit with the time-dependent 
signal 

Wt) = *o/2 + $ i + Trgl($ / ,r,t), (I) 

where is the initial flux detuning and Trgl($/,r, t) 
is the triangle signal parameterized by $/ and r, which 
correspond to the final value of the flux sweep and the 
time width respectively(shown in Fig.flSk)). The explicit 
expression is 



Trgl($/, r, t) 



kt (0 < t < i) 

k(t-T) (L<t<T) ' [Z) 

where k is the sweep rate 

k = -. (3) 

T 

First we consider the two lowest states |L0) and \R0), 
that is, to treat the qubit as a pure two-level system. 
The reduced Hamiltonian is 
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(4) 



±Q(& ext ) are energy frequencies of the ground states in 
two wells; A is the tunneling frequency between the two 
states. In following text we set H = 1. When we add 
signal to § e xt, the Hamiltonian is time-dependent. We 
use density matrix to calculate the populations in the 
states. The density matrix of the qubit is 



W n (t) W 12 (t) 
W 21 (t) W 22 (t) 



(5) 



Wu and W 22 are populations in |i0) and \R0) respec- 
tively, and obey the unity condition W\\+W 22 = 1. W\ 2 
and W 2 i, which are complex conjugate, mark the coher- 
ence of the two states. The time-evolution of the density 
matrix satisfies the celebrated Louville equation 

■ dp(t) 



dt 



-iTp(t) + [H red J(t) 



(6) 



The decaying rate tensor f will be ignored because the 
signal width is much shorter than the decaying time. 

We initialize the qubit on \L0), which means that the 
initial condition is 



1 




(7) 



Solve the Louville equation Eq.© to obtain the value of 
Wn(t) at t = t, the population in \LQ) when the sweep 
signal ends. Choosing ranges for the signal parameters 
$/ and r, and solving the Louville equation for all com- 
binations of signal parameters within their ranges, we 
obtain a map showing the final Wu versus the two pa- 
rameters. 
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III. RESULTS AND ANALYSIS 

Shown in Fig. ([2]) are our simulations. In this case, we 
preset the spectrum with an energy slope of 2 GHz/m$ 
and A of 2 GHz. For the signal we have $i = —5 m$ 0l T 
ranging from 0.01 to 4 ns and $ / from -2 to 10 m<f>o- The 
final population in the initial state \L0) oscillates along 
both axes, demonstrating the constructive interference 
and the destructive interference. The oscillating periods 
are determined by the energy parameters of the qubit, so 
we can extract the qubit spectrum from the interference 
pattern. 




It is easy to obtain an analytical solution, which can 
quantitatively elaborate the result and give specific steps 
to extract data. Sweeping through the anticrossing incurs 
Laudau-Zener transitions, giving rise to populations in 
l-RO). When sweeping back to the anticrossing, a phase 
difference of the two states \LQ) and \R0) is accumulated: 



9 = 



{vi(t) - vo(t))dt 



(8) 



Pq (f ) and v\ (t) are given by the eigenvalues of Hamilto- 
nian Q: 



(9) 



r* is the effective width of the signal, equating the inter- 
val of the two times of passing the anticrossing (shown 
in Fig. (SO): 



2$ 



(10) 



k <&/ - $j 

Notice that 17(f) varies linearly with time 

17(f) = J(*«rt-$o/2) = Z(fcf-$ /2)(0 < f < r/2), (11) 

where I is the slope of the spectrum. Substituting it into 
the integral Eq.([S]) and resetting the time parameter, we 
have 



o 



(vx(t) - vo{t))dt 



r*/2 



(12) 



2 / 2y / (lkt) 2 + A 2 dt 
la 



= Ikr* 



, A N 2 / T* , ; 



A 2 , , Ikr* 
-2 — In — - 
Ik 1 2A 



1 



,lkr* 



(13) 



Using Eq.Q and Eq. (fl0|) to replace k and t* with the 
original signal parameters, we can write the above equa- 
tion as 
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1$, 



■In 



(14) 




Amplitude, A(m<I) ) 




Amplitude, A (m* ) 



FIG. 2: (color online) (a)Illustration of the qubit spectrum. 
Only the two lowest levels |L0) and |J?0) are taken into ac- 
count. Preset parameters: the slope is 2 GHz/m^o and A is 2 
GHz. The inset shows the details of the signal, with the time 
width r and the final value $/. "3>i is the initial flux detuning 
and set to be -5 m$o in this simulation, r* is the interval of 
two times of passing through the anticrossing. (b)The sim- 
ulation result of LZS interference. The qubit is driven with 
the signal illustrated in (a). The population in |L0) begins to 
oscillate when the anticrossing is reached. The data at a large 
amplitude, e.g. along the white dashed line, can be used to 
calculate the slope of the energy spectrum. (c)The data on 
the white dashed line: population in |L0) versus the signal 
width r when the amplitude $/ is fixed at 8 m$o. The pe- 
riod of the oscillation is about 0.7ns. As a fitting parameter, 
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It is this phase difference that causes the interference. It 
can proved that the population in the initial state oscil- 
lates sinusoidally with 

In the limit of large amplitude driving we have 



1*1 
A 



»1, 



thus 



w « — r 



(15) 



from which I can be extracted as the fitting parameter. 

To test our method, we extract a section in Fig. (j^t 1 ) 
along the white dashed line. This section (see Fig.([5J:)) 
shows the relation between the population in |L0) and the 
signal width r at a fixed amplitude that is large enough 
to validate the approximation in Eq. (1151) . The oscillation 
period T is about 0.7 ns. According to Eq.(fT5"j). 



2tt 



2tt 
07 



ns 



Note that <&/ is fixed at 8 m$ n and 



-5 m<I>o, the 

slope is 1.83 GHz/mc&o- Comparing with our preset pa- 
rameter I — 2 GHz, we reach a precision of at least 90%. 
In the extremely amplitude regime, Eq. (|15j) can be again 
approximated as 



l<t> 



P 



Thus, 



2tt 



./'■ 



(16) 



(17) 



To verify this prediction, we extend the calculation range 
of and make a discrete one dimensional Fourier 
transformation to the interference pattern. The one di- 
mensional Fourier transformation diagram is shown in 
Fig.([2ji). We can clearly observe the linear relation at 
extremely large values of $/ . Actually, when $/ is larger 
than 8 m$o, the higher levels should be taken into ac- 
count, as is the case discussed in the following section. 

After evaluating Z, we can make a further step to ex- 
tract A from Eq. (fT4")) . We only need to choose two points 
in the lower region of Fig.([2jD), for instance: 
$ f = 1.0 m$ , r = 3.85 ns, Population = 0.00; 
$/ = 2.0 m$ ,T = 3.85 ns, Population = 0.93. 
The population on |L0) is equal to22 

1 



-(1 + cosp), 



(18) 



in which ip can be expressed in terms of A as in Eq. (|14p . 
So we plot Eq. ([T5| with A as the x-axis and find the 
value of A which fits the both points we have chosen. 
The result is shown in Fig ([3]). The only coincident result 
is A = 2 GHz, which is also precisely in accordance with 
our preset parameter. 




FIG. 3: Plots of the function in Ea. (fl8|) with two sets of 
parameters. We find that at A = 2 GHz, both values fit 
well with the simulation, indicated with the arrow. Thus we 
obtain the spectrum parameter A. 
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FIG. 4: (color online) Illustration of the spectrum of the three 
levels in question. A12 and A13 mark the sizes of the energy 
gaps at the anticrossings. The bold arrow illustrates the large 
amplitudes that sweep through both anticrossings. 



IV. MULTI-LEVEL SYSTEM 

If we drive the qubit with a signal whose amplitude 
is large enough to reach another anticrossing in higher 
energy levels, the interference pattern emerges with in- 
teresting characteristics, especially for some sets of pa- 
rameters. The simulation approach is quite similar to 
the case in Section Hill We consider the three states \L0), 
\RQ) and \R1), as shown in Fig.Q. Under the basis com- 
posed of these three states, the reduced Hamiltonian can 
be written as 



ui A12 Ai 3 

Hred = H | A 12 U) 2 

Ai 3 uj 3 



(19) 



Still we will set H = 1. 0J1, oj?, and W3 are energy frequen- 
cies corresponding to the three states in question. A12 
and A13 are tunneling frequencies between |L0) and the 
other two states respectively. They mark the scales of 
the energy gaps at the two anticrossings. \R0) and \R1) 
are not correlated. 
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Now we have a 3 x 3 density matrix 

fW u W 12 W 13 \ 
p = W 21 W 22 W 23 . (20) 

V w 31 w 32 w 33 J 

The calculation procedures will be the same as in Sec- 
tion |nH We have changed the ratio A13 : A12 and simu- 
lated the interference in three cases, each of them having 
unique characteristics. Now we explain the three cases 
in details. 

Case 1 (Fig.([SJi)): A 13 = 10 GHz, A 12 = 1 GHz. 
We have the ratio Ai 3 : A i2 as large as 10 and ob- 
tain an interference pattern that looks similar to the one- 
anticrossing case. The reason is that the first anticross- 
ing, marked by A12, is too small to affect the results. 
There is no apparent oscillations when < $/ < 8 m$o, 
since strong transitions keep occurring at the vicinities as 
if there were no energy gaps. So, the interference fringes 
are caused by the bigger anticrossing, whose location can 
be easily inferred from the left edge of the first fringe. 
In this case we find it at 8 m$o, in agreement with the 
spectrum in Fig.Q. Moreover, the parameters related 
to this anticrossing can also be extracted, following the 
steps introduced in Section lHll However, the information 
of the smaller anticrossing cannot be obtained because its 
effects are not revealed in this case. 
Case 2 (Fig.([5p)): A 13 = 8 GHz, A 12 = 2 GHz. Then 
we lower the ratio A13 : A12 to 4 and then observe a 
significant difference. The stripes are distorted and ir- 
regular when the second anticrossing is reached. The 
distortion gets slighter as the signal width increases and 
the pattern becomes close to that in the first case. This 
can be anticipated because increasing width means slower 
sweeping and thus gives rise to weaker transitions at the 
second anticrossing. In the lower part of the map, the 
pattern resembles that in the one-anticrossing situation, 
since all that contributes is mainly the first one. Here we 
can obtain overall information of the energy spectrum, 
because the two anticrossings are of commensurate sizes 
that dual effects are observed. We will particularly dis- 
cuss this case later since it is of special interests. 
Case 3 (Fig.©)): A 13 = 2 GHz, A 12 = 8 GHz. Now 
we set the ratio A i3 : A i2 = 1/4. We encounter a result 
similar to a one-anticrossing case again. Like the first 
case, the reason is that the contribution of the smaller 
anticrossing is rather trivial. It is the first bigger an- 
ticrossing that primarily contributes in the small-width 
region. In the figure we find clear interference fringes 
from = 0, where the first anticrossing is located. And 
the details concerning this anticrossing can be obtained 
in the same way as in Section IIIII 

In all the three figures there are very delicate ripples, 
which result from the high precision of our calculation. 
While in experiments, the ripples can hardly observed 
due to noises, decaying or measurement limits. 

Among the three cases we have discussed, the second is 
of special interests and significance because it is the clos- 
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FIG. 5: (color online) Simulation results of LZS interference 
with two anticrossings involved, marked by A13 and A12 re- 
spectively. (a)In this case, we have A13 : A12 = 10. The inter- 
ference fringes primarily result from the bigger anticrossing, 
which can be located by the left edge of the first fringe. Its 
parameters can be also obtained by analyzing the pattern as a 
one-anticrossing case, (b) Commensurate energy gaps are set. 
The ratio A13 : A12 is 4. The map can be divided into three 
regions by two characteristic sweep rates fci2 and fci3, shown 
with dot lines. In Region(3), distorted interference fringes are 
observed, due to the cooperative effects of both anticrossings. 
While in Regions(l) and (2), the pattern is similar to the one- 
anticrossing case, since one of the anticrossings loses its effect. 
Thus the information of the anticrossings can be inferred from 
these two regions respectively. (c)Here A13 : A12 = 1:4, re- 
sulting in a pattern looking like (a). The big anticrossing 
overwhelms the small one in most area. Its information can 
be obtained as well. 
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est model to the real situation. This is a typical multi- 
level system with two anticrossings, on which our LZS 
spectroscopy method is applicable as well. First of all, 
the locations of the anticrossings can be inferred from the 
pattern. Similarly, the left edge of the first fringe tells 
the location of the first anticrossing(at m$o), and the 
beginning of the distortion marks the second anticrossing 
(at 8 m$o)- For further elaborations, we divide the map 
into three regions with two dashed lines representing two 
characteristic sweep rates ku, such that 



referring to Landau-Zener transition rate formul a 31 ' 33 . 
In Regional ) where k < ki2, the pattern is mainly pro- 
duced by the first anticrossing A12, so we can calculate 
the energy slope and the value of A12. Similarly, in 
Region(2) where k > fci3, the information of the sec- 
ond anticrossing A13 can be obtained. Region(3) clearly 
demonstrates the cooperative effects of both anticross- 
ings and the transition from Region(l), dominated by 
A12, to Region(2), dominated by A13. Along the direc- 
tion of increasing sweep rates, the distortion results in 
denser stripes. It is because that the bigger anticrossing 
contributes more when sweeping is faster, so the accumu- 
lated phase changes more rapidly, according to Eq. fTB")) . 

There have been some experimental results that qual- 
itatively prove our simulations. Recently, some re- 
searchers have reported that the spectrum of an rf 
SQUID is modified by some spurious oscillators. One of 
the significant works was carried out by Simmonds et al. 
in 20042i. The spurious oscillators, later interpreted as 
two-level system (TLS), essentially and universally exist 
in Josephson devices. They transform the energy struc- 
ture of a phase qubit and lead to multiple anticrossings 
by coupling with the qubit. One of the latest works 
on such configuration is held by Sun et al£^, who per- 
formed Landau-Zener interferometry in a TLS-coupled 
phase qubit. The spectrum of the qubit shows two TLSs 
coupled to the system, leading to two anticrossings on the 
left side of § e xt = 3>o/2- This two-anticrossing structure 
is similar to the one we have discussed, although their ori- 
gins are different. Sun et al. used triangle pulses to drive 
the system, and observed the interference fringes charac- 
terized by the distortion, which has also been revealed in 
the second case of our simulations. Notably, Shevchenko 
et al. also have made some deep investigations into this 
issue very recently 3 ^. 

Additionally, we want to point out that the method of 
LZS spectroscopy suggests a possibility of coherent ma- 
nipulations on a flux qubit. Since the whole process has 
a time scale of a very short triangle pulse, it ensures that 
the operation is performed within the decoherence time 
of the qubit. Moreover, with different sweep rates, we 
can switch the two anticrossings on and off respectively. 
The switching sweep rates are exactly the two we have 
used to divide the map. We can interpret the three re- 
gions from an alternative view: 



(1) fc < k\2- Transitions occur mainly at the first anti- 
crossing but are ignorable at the second, as if the second 
anticrossing is switched off. 

(2) fc > fci3. Now the first anticrossing is switched off 
while the second is on, because transitions occurring at 
the first anticrossing are very strong but moderate at the 
second. 

(3) fci2 < k < /C13. This situation corresponds to the 
region where the two anticrossings are both switched on. 

Therefore, the anticrossings act like tunable beam 
splitters with transmission coefficients changing continu- 
ously from null to unit 3 ^. By adjusting the sweeping rate 
carefully, we can control the population in the excited 
states, supplying an alternative method to manipulate 
the qubit. 



V. CONCLUSIONS 

We propose a new method of measuring the energy 
spectrum of a superconducting flux qubit. By sweeping 
the qubit through the energy anticrossings with a linear 
triangle pulse, we obtain the LZS interference patterns. 
Then by fitting the patterns with analytical equations 
we can extract the information of the energy spectrum, 
including the slope of the spectrum and the magnitudes 
of energy gaps at anticrossings. We have demonstrated 
this method with numerical simulations. It can be a 
more convenient and efficient method with a precision 
over 90%, especially in measuring a spectrum with high 
energies, in which conventional frequency spectroscopy 
has difficulties. 

Another area where our method may be useful is in 
the adiabatic quantum computation 3 ?* 3 ^, in which mi- 
crowaves are totally unnecessary. If we use this linear 
spectroscopy to calibrate the qubits, the measurement 
setup may be simplified and better shielded from the ex- 
tra noise, which is otherwise introduced from microwave 
lines. 

This method also can be applied to other quantum 
systems in which rf field coupling may be too weak to 
generate detectable population transitions, so we can try 
linear signals instead of rf fields to perform spectroscopy 
and manipulations. 

In addition, since the operation is done within the de- 
coherence time of a qubit, it can be utilized to realize 
coherent manipulations of the qubit. We have also dis- 
cussed the approach of controlling the functions of the 
anticrossings. For our simulations qualitatively agree 
with related experimental results, we think that our work 
would stimulate more investigations in this field. Mean- 
while, this approach can be extended to multi-qubit sys- 
tems, whose energy spectra are of analogous structures. 
Therefore, the approach might also shed light on large- 
scale controllable quantum computation in future. 
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